Casimir interaction between plane and spherical metallic surfaces 
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We give an exact series expansion of the Casimir force between plane and spherical metallic 
surfaces in the non trivial situation where the sphere radius R, the plane-sphere distance L and 
the plasma wavelength Ap have arbitrary relative values. We then present numerical evaluation 
of this expansion for not too small values of L/R. For metallic nanospheres where R,L and Ap 
have comparable values, we interpret our results in terms of a correlation between the effects of 
geometry beyond the proximity force approximation (PFA) and of finite reflectivity due to material 
properties. We also discuss the interest of our results for the current Casimir experiments performed 
with spheres of large radius R^ L. 



The Casimir force is a striking macroscopic effect of 
quantum vacuum fluctuations which has been seen in a 
number of dedicated experiments in the last decade (see 
for example and references therein). One aim of the 
Casimir force experiments is to investigate the presence of 
hypothetical weak forces predicted by unification models 
through a careful comparison of the measurements with 
quantum electrodynamics predictions. This aim can only 
be reached if theoretical computations are able to take 
into account a realistic and reliable modeling of the ex- 
perimental conditions. Among the effects to be taken 
into account are the material properties and the surface 
geometry, these effects being also able to produce phe- 
nomena of interest in nanosystems [s!, f3| . 

A number of Casimir measurements have been per- 
formed with gold-covered plane and spherical surfaces 
separated by distances L of the order of the plasma wave- 
length (Ap ~ 136nm for gold), making material proper- 
ties important in their analysis 5]. As those measure- 
ments use spheres with a radius R ^ L, they are com- 
monly analyzed through the Proximity Force Approxi- 
mation (PFA) [6] , which amounts to a trivial integration 
over the sphere-plate distances. An exception is the Pur- 
due experiment dedicated to the investigation of the ac- 
curacy of PFA in the sphere-plate geometry [Ti] , the result 
of which will be given as a precise statement below. 

In the present letter, we give for the first time an ex- 
act series expansion of the Casimir force between a plane 
and a sphere in electromagnetic vacuum, taking into ac- 
count the material properties via the plasma model (see 
Fig. 1). We present numerical evaluation of this expan- 
sion which are limited to not too small values of L/R, 
because of the multipolar nature of the series. We show 
below that these new results lead to a striking correlation 
between the effects of geometry and imperfect reflection 
when evaluated for nanospheres, with R, L and Ap hav- 
ing comparable values. In the end of this letter, we also 
discuss the interest of these results for the Casimir ex- 
periments performed with large spheres i? ^ L Q ■ 

Our starting point is a general scattering formula for 
the Casimir energy 0. Using suitable plane- wave and 
multipole bases, we deduce the Casimir energy £ps be- 




FIG. 1: The geometry of a sphere of radius R and a flat plate 
at a distance L (center-to- plate distance £ = L -\- R); both 
mirrors are covered with a metal characterized by a plasma 
wavelength Ap. 



tween a plane and a spherical metallic surface in elec- 
tromagnetic vacuum. The multipole series expansion is 
written in terms of Fresnel reflection amplitudes for the 
plate and Mie coefflcients for the sphere, and it is valid 
for arbitrary relative values of the sphere radius R, the 
sphere-plate distance L and the plasma wavelength Ap. 
For the sake of comparison with experiments, we assume 
Ap ~ 136nm for both, the sphere and the plate. We 
occasionally also consider the limit Ap — > 0, where the 
formula reduces to the case of perfect reflectors in elec- 
tromametic vacuum, for which results were obtained ve- 
cently 

In the following, we discuss the force JTpg = —d£ps/dL 
as well as the force gradient t/pg = —dJ-ps/dL which was 
measured in the experiment [7j. We write the results 
deduced from the scattering formula as products of PFA 
estimates by beyond-PFA correction factors pp and po- 
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The PFA estimates J-pg^ and Qpg^ are proportional re- 
spectively to the energy and force calculated between two 
planes. They are written as products of ideal Casimir ex- 
pressions and factors ?7e and r]p accounting for the effect 
of imperfect reflection Q . 
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The beyond-PFA correction factors pp and pc appear- 
ing in ([T]) are the important quantities for what foUows. 
For experiments performed with large spheres of radius 
L, the deviation from PFA is small (pp — !)■ Even 
in this limit, it remains important to specify the accu- 
racy of PFA in order to master the quality of theory- 
experiment comparison [To^. This can be done by intro- 
ducing a Taylor expansion of the correction factors at 
small values oi L/R 
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The only experimental result available on this topic 
may be stated as a bound on the Pc factor, namely 
|/3g| < 0.4. On the theoretical side, analytical as well as 
numerical calculations of this slope have been obtained 
for scalar field models [H [H, [M [H, [3 . For the situ- 
ation met in experiments, with a plane and a sphere in 
electromagnetic vacuum, an estimation technique has re- 
cently been proposed where the slope is deduced from a 
polynomial fit of the numerical values obtained at inter- 
mediate values of L/R [1, [13]. The slope obtained in this 
manner is much larger (^8 times larger) than expected 
from scalar field models [lO] . As a consequence, the value 
of Pg falls out of the bound of 0], in contrast with the 
scalar prediction which lies within the bound. More pre- 
cise statements on this point will be given below. 

On the other hand all these results correspond to per- 
fect reflection, whereas the experiment Q was performed 
with gold-covered surfaces. The apparent contradiction 
noticed in the preceding paragraph may thus be cured if 
the value of Pg differs for metallic and perfect mirrors, 
that is also if the effects of geometry and finite reflectiv- 
ity are correlated. We show in the sequel of the letter 
that this is indeed the case. 

We start from the formula for the Casimir energy fpg 
between two scatterers in vacuum Isl 
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In the geometry depicted on Fig. 1 with a sphere of radius 
R, a plate, and a sphere-plate separation L along the z- 
axis (center-to-plate distance £ = L + R), TZs and TZp 
represent the reflection operators for the spherical and 
the plane scatterers, respectively. They are evaluated 
with reference points placed at the sphere center and at 
its projection on the plane, respectively. The operator 
g-icc describes the one-way propagation between these 
two reference points. ^ is the imaginary fleld frequency 
integrated over the upper imaginary axis. 

In order to evaluate explicitly this expression, we use 
two mode decompositions. The first one is a plane-wave 
basis |k, with k the transverse wavevector paral- 

lel to the xy plane, p = TE, TM the polarization, and 
(j) = ±1 for rightward/leftward propagation directions. 
It is well adapted to the description of free propagation 



and reflection on the plane: the propagation operator 
g~K.c jg diagonal with matrix elements e'^^ such that 
K ~ y^|^/(?+fc2 (^j^ = while reflection on the plane 
preserves all plane-wave quantum numbers but (j). The 
non zero elements of TZp are the standard Fresnel reflec- 
tion amplitudes r^. Given values of k(fc, ip) and (j) = ±1 
define a direction in reciprocal space corresponding to 
the azimuthal angle ip and a complex angle 9^ such that 
sin 6'='= = -i^ and cos^='= = 

The second basis, which is adapted to the spherical 
symmetry of TZs, is a multipole basis \£mP)^, with £{£+l) 
and m the angular momentum eigenvalues {i — 1,2, 
m — and P = E,M for the electric and mag- 

netic multipoles. By rotational symmetry around the 
2;-axis, A4 commutes with Jz. Hence it is block diagonal, 
with each block corresponding to a common value 

of m and yielding a contribution fp™^ to the Casimir 
energy Spg (opposite values ±to provide identical contri- 
butions). The contribution i?p™'' is written as in ^ with 
Ai replaced by the block matrix 
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Each block in this matrix is the sum of TE and TM 
contributions M(™)(Pi, Pa) = J2p a4"^(Pi, P2). The di- 
agonal blocks are written as 
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ai{i£/) and 6^(1^) are the Mie coefficients [T^l for electric 
and magnetic multipoles. A and B are matrices which 
do not depend on the radius nor on the refractive index 
of the sphere and are written in terms of the spherical 
harmonics Yt^rn{0,^p = 0) and the finite rotation matrix 
elements df^^^,{e) = {£,m\e~^^Jy\£,m') [l8| 
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Similar expressions are found for the nondiagonal blocks, 
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with the matrices A and B replaced respectively by 

xY,,„,{d-) rp{k)e-^''^ (7) 

In order to go further, we assume the materials to 
have a dielectric response described by the plasma model 
e(iO = 1 + '^p/C^j with Up the plasma frequency and 
Ap = 2ttc/ujp the plasma wavelength. Although the for- 
malism easily allows for different values of Ap for both 
surfaces, we take a common value as in the recent ex- 
periment We calculate the Casimir energy £ps and 
deduce the force ^ps and gradient ^ps, both quantities 
being functions of the 3 length scales R, L and Ap. The 
case of perfect reflection fld\ can be recovered as the limit 
Xp R,L (see ^9|] for the opposite non retarded limit). 
A large distance limit may also be taken as Ap, R ^ L. 
Its result reduces to the Rayleigh expression (2(| in the 
case {R < Ap) or to 3/2 of it [1, [lO| in the case (Ap ^ R). 

As already discussed, the PFA expression is also con- 
tained in our general result, and it is recovered asymptot- 
ically for i? L. In the following, we discuss the results 
of numerical computations of the ratios pf,g defined in 
([1]) which measure the deviation from PFA. For dimen- 
sionality reasons pf,g are functions of two dimensionless 
parameters built upon L, R and Ap (t?e,f are functions 
of L/Xp only f^) and they approach unity at the PFA 
limit L/R ^ 1. Their numerical computation is done af- 
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FIG. 2: Upper graph : variation of pp as a function of L/R, 
for a nanosphere of radius R = lOOnm; the solid green line 
corresponds to gold-covered plates (Ap = 136nm) and the 
dashed red line to perfect reflectors. Lower graph : variation 
of pg as a function of L/R, with the same conventions as on 
upper graph. The decreases at low values of L/R represent 
a numerical inaccuracy due to the limited value of ^max = 24 
[Colors online]. 



ter truncating the vector space at some maximum value 
^max of the orbital number i. As a consequence of the 'lo- 
calization principle' [llj, the results are accurate only for 
R/L smaller than some value which increases with £niax- 
At the moment, our numerical calculations are limited to 
^max = 24, allowing us to obtain accurate results down 
to L/R ~ 0.2 but not in close vicinity of the PFA limit. 

This method gives new and interesting results, in par- 
ticular for nanospheres having a radius R with the same 
order of magnitude as the plasma wavelength Ap. In this 
case, we can perform accurate calculations for L hav- 
ing a comparable magnitude, and thus explore the rich 
functional dependence of Pf,g versus two dimensionless 
parameters built up on L,R and Ap. Fig. [2] shows the 
results obtained for pp and pQ with metallic and per- 
fect mirrors. Clearly the deviation from PFA calculated 
for metallic mirrors differs markedly from that already 
known for perfect mirrors. For small values of L/R the 
violation of PFA for the Casimir force and gradient turns 
out to be less pronounced for metallic mirrors than for 
perfect mirrors, while for large values of L/i? it is more 
pronounced. 

However, at values L/R ~ 0.2 we find a clear correla- 
tion between geometry and finite reflectivity effects, mak- 
ing therefore measurements with nanospheres at small 
plate-sphere separations particularly interesting. This 
non trivial interplay becomes evident when a polynomial 
fit of the numerical values of pf,g is used for inferring 
the behaviour at smaU values of L/R On Fig. [3] 

we plot the quartic polynomial fits of the function pc for 
the two cases of gold-covered and perfect mirrors. The 
curves were obtained by finding the best-fit of the nu- 
merically computed values of pc (crosses on Fig. [3]) in 
the window 0.4 < L/R < 0.8 (circled crosses on Fig. [3]) 
in the set of quartic polynomials (Taylor expansion de- 
fined as in (12) and truncated at fourth order). The left- 
hand bound of the window is fixed by the requirement of 
using only points accurately calculated with fmax = 24 
while the righthand bound is determined by the trunca- 
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FIG. 3: Quartic polynomial fit of the function pg{L/R), for 
a nanosphere of radius R = lOOnm; the solid green line cor- 
responds to gold-covered plates and the dashed red line to 
perfect reflectors. The crosses represent numerically evalu- 
ated points and the circles indicate those points which are 
used for the fit [Colors online]. 
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tion at fourth order of the Taylor expansion. The best-fits 
correspond to the following polynomials for gold-covered 
(GM) and perfect (PM) mirrors respectively {x = L/R) 

GM : 1 - 0.2072: - O.SSOx^ + OMbx^ - 0.249a;^ 
PM : 1 - 0.4832: 4- 0.297a;2 - 0.221x3 + O.OSOa:^ 

The two fits are clearly different and this in particular the 
case for the values obtained for the slope at L/R — 0. 
The slope (/3g ~ —0.21) obtained for gold mirrors dif- 
fers by more than a factor 2 from the one (/3g ~ —0.48) 
obtained for perfect mirrors. This is related to the bend- 
ing of the curve for gold mirrors at small L/R, which 
describes the effect of imperfect reflection in the beyond- 
PFA factor pQ and has to be contrasted with the un- 
bent curve for perfect mirrors. For the same reason, we 
observe that the slope obtained for gold mirrors is less 
stable under the variation of the conditions of the best- 
fit procedure than that for perfect mirrors. To appre- 
ciate the meaning of the bending let us recall that the 
slope obtained for perfect mirrors in an electromagnetic 
vacuum is ~8 times larger than expected from scalar 
computations [H, [lB| and one cannot but notice that 
it lies outside the bound \/3g\ < 0.4 of Q- In con- 
trast, the slope obtained for metallic mirrors lies within 
the bound. Let us emphasize that there is no contra- 
diction between the results presented here (obtained for 
nanospheres with R = lOOnm) and the experiments (per- 
formed with spheres having i? > a few tenths of /im) . 

For spheres with large radii {L/R > 0.2) the beyond- 
PFA factors pf,g have the same values for gold-covered 
and perfect mirrors, because the value of L is much larger 
than Ap. If we extracted a slope from these results, we 
would obtain a value close to that of perfect mirrors, thus 
lying outside the bound of Q. However, the arguments 
discussed before show that one should refrain from doing 
so. Indeed, a bending of the curve has to be expected in 



this case too, for values of L becoming comparable to Ap 
and thus much smaller than R. In contrast, this bending 
has no reason to appear for perfect mirrors since there 
is no length scale like Ap in this case. If the bending is 
similar for large and small spheres, it may turn out that 
the slope for gold-covered mirrors meets the bound [3| 
while that for perfect mirrors does not. 

To sum up our results, we have written a new and ex- 
act expansion for the Casimir force between plane and 
spherical metallic surfaces in electromagnetic vacuum. 
The results go beyond the proximity force approximation, 
and show a clear correlation between the plane-sphere 
geometry and the material properties of the metallic sur- 
faces. They constitute a new step in the direction of 
accurate comparisons between Casimir experiments and 
QED theoretical predictions. More work is needed to ob- 
tain exact results for the Casimir force between a metal- 
lic sphere and plate in the so far experimentally explored 
parameter region of L/R ~ 0.01, using for example dif- 
ferent approaches based on semiclassical methods. Our 
results also indicate a complementary way to observe de- 
viations from PFA and the interplay between geometri- 
cal and reflectivity effects in new experiments performed 
with nanospheres. 
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